Theory of interacting electrons on the honeycomb lattice 
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The general low-energy theory of electrons interacting via repulsive short-range interactions on 
graphene's honeycomb lattice at half filling is presented. The exact symmetry of the Lagrangian 
with local quartic terms for the Dirac four-component field dictated by the lattice is identified as 
D2 X (/c(l)xtime reversal, where D2 is the dihedral group, and r/c(l) is a subgroup of the SUc{2) 
"chiral" group of the non- interacting Lagrangian, that represents translations in Dirac language. 
The Lagrangian describing spinless particles respecting this symmetry is parameterized by six inde- 
pendent coupling constants. We show how first imposing the rotational, then Lorentz, and finally 
chiral symmetry to the quartic terms - in conjunction with the Fierz transformations - eventually 
reduces the set of couplings to just two, in the "maximally symmetric" local interacting theory. We 
identify the two critical points in such a Lorentz and chirally symmetric theory as describing metal- 
insulator transitions into the states with either time-reversal or chiral symmetry being broken. The 
latter is proposed to govern the continuous transition in both the Thirring and Nambu-Jona-Lasinio 
models in 2+1 dimensions and with a single Dirac field. In the site-localized, "atomic", limit of the 
interacting Hamiltonian, under the assumption of emergent Lorentz invariance, the low-energy the- 
ory describes the continuous transitions into the insulator with either a finite Haldane's (circulating 
currents) or Semenoff 's (staggered density) masses, both in the universality class of the Gross-Neveu 
model. The simple picture of the metal-insulator transition on a honeycomb lattice emerges at which 
the residue of the quasiparticle pole at the metallic, and the mass-gap in the insulating phase both 
vanish continuously as the critical point is approached. In contrast to these two critical quantities, 
we argue that the Fermi velocity is non-critical as a consequence of the dynamical exponent being 
fixed to unity by the emergent Lorentz invariance near criticality. Possible effects of the long-range 
Coulomb interaction, and the critical behavior of the specific heat and conductivity are discussed. 



I. INTRODUCTION 

Two-dimensional honeycomb lattice of carbon atoms 
may be viewed as the mother of all other forms of car- 
bon. Its crucial electronic property, which arises as a 
consequence of the absence of the inversion symmetry 
around the lattice site, is that the usual Fermi surface is 
reduced to just two points. The electronic dispersion may 
be linearized around these two points, after which it be- 
comes isotropic and dependent on the single dimension- 
ful parameter, Fermi velocity vp « c/300. The pseudo- 
relativistic nature of the electronic motion in graphene 
has since its synthesis placed this material at the center 
stage of condensed matter physics. Many qualitatively 
novel phenomena that take, or may take place in such a 
system of " Dirac" electrons are actively discussed in the 
rapidly growing literature on the subjectii 

In this paper we discuss the low-energy theory and the 
metal-insulator quantum phase transitions of the inter- 
acting Dirac electrons on the honeycomb lattice, build- 
ing upon and expanding significantly the earlier work 
by one of uSfS In the first approximation, all weak in- 
teractions of Dirac electrons in graphene may be ne- 
glected at half filling, when the Fermi surface consists 
of the Dirac points. This is because short-range inter- 
actions are represented by local terms which are quartic 
in the electron fields, which makes them irrelevant near 
the non-interacting fixed point by power counting. The 
same conclusion turns out to apply to the long-range tail 
of the Coulomb interaction, which remains unscreened 
in graphene, although only marginally soi^>^ Neverthe- 



less, if strong enough, the same interactions would turn 
graphene into a gapped Mott insulator. As an example, 
at a strong on-site repulsion the system is likely to be 
the usual Neel antiferromagneti^ It is not a priori clear 
on which side of this metal-insulator transition graphene 
should be. With the standard estimate for the nearest- 
neighbor hopping in graphene oi t — 2.5eV and the Hub- 
bard interaction oi U « 7 — 12eV, it seems that the 
system is below yet not too far from the critical point es- 
timated to be at U/t « 4— Si^iii^ If sufficiently weak, the 
electron-electron interactions only provide corrections to 
scaling of various quantities, which ultimately vanish at 
low temperatures or frequencies. At, what is probably a 
more realistic, an intermediate strength, the flow of in- 
teractions and the concomitant low-energy behavior may 
be influenced by the existence of metal-insulator criti- 
cal points. It is possible that some of the consequences 
of such interaction-dominated physics have already been 
observed in the quantization of the Hall conductance at 
filling factors zero and onej^ii^iiiii^ii^iii^ As we argued 
elsewhere, the anomalously large value of the minimal 
conductivity in graphene^^ may be yet another conse- 
quence of the Coulomb repulsion between electrons J^iii 

The above discussion raises some basic questions. 
What is the minimal description of interacting electrons 
in graphene at " low" energies? What is the symmetry of 
the continuum interacting theory, and how does it con- 
strain the number of coupling constants? What kinds of 
order may be expected at strong coupling, and what is 
the nature of the metal-insulator quantum phase transi- 
tion? In this paper we address these and related issues. 
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FIG. 1: Two axis of symmetry of the low-energy theory of 
graphene, in real space. The Dirac points in this coordinate 
frame are at ikK = (1, 0)(47r/3a), i. e. along the A-axis. 

In the rest of the introduction we give a preview of our 
main results. 

The simplest prototypical system that exhibits the 
physics of interacting Dirac fermions which we seek to 
understand is the collection of spinless electrons inter- 
acting via short-range interactions, at half-filling. For 
present purposes an interaction may be considered as 
"short-ranged" if its Fourier transform at the vanishing 
wavevector is finiteJ^ The least irrelevant quartic terms 
one can add to the non-interacting Dirac Lagrangian will 
then be local in space-time, and of course quartic in terms 
of the four-component Dirac fields that describe the elec- 
tronic modes near the two inequivalent Dirac points at 
wavevectors ztK at the edges of the Brillouin zone. The 
most general local quartic term in the Lagrangian would 
be of the form 

L^nt = (*^(x,r)A/i*(x,T))(*t(f,r)M2*(f,r)), (1) 

where Mi and M2 are four-dimensional Hermitian matri- 
ces. The symmetry alone, however, immediately drasti- 
cally reduces the number of independent couplings from 
the apparent 136 to just fifteen. Although the point 
group of the honeycomb lattice is Cg^, the exact spatial 
discrete symmetry of the Lagrangian is only the dihedral 
group D2, or the vierergruppe, which consists of the re- 
flections through the two coordinate axis shown in Fig. 
1, and the inversion through the origin. Such a small 
symmetry results from the very choice of two inequiva- 
lent Dirac points out of six corners of the Brillouin zone, 
which reduces the symmetry to the simple exchange of 
the two sublattices (reflection around A axis), the ex- 
change of Dirac points (reflection around B axis), and 
their product (the inversion through the origin). D2, the 
time-reversal, and the translational invariance are shown 
to leave fifteen possible different local quartic terms in 
the Lagrangian. 

Fortunately, not all of these still numerous quartic 
terms are independent, and there are linear constraints 
between them implied by the algebraic Fierz identities.J^ 
The Fierz transformations are rewritings of a given quar- 
tic term in terms of others, and we provide the general 
formalism for determining the number and the type of 



independent quartic couplings of a given symmetry. For 
the case at hand we find that spinless electrons interact- 
ing with short-range interactions on honeycomb lattice 
are in fact described by only six independent local quar- 
tic terms. The inclusion of electron spin would double 
this number to twelve. 

The linearized noninteracting Lagrangian for Dirac 
electrons, 

Lo^^{x,T)j^df,^{x,T) (2) 

as well-known, exhibits the Lorentz and the global 
SUci^) ("chiral") symmetry. The latter, generated by 
{73, 75, 735}, with 735 = -«7375, is nothing but the "rota- 
tion" of the " pseudospin" , or " valley" , corresponding to 
two inequivalent Dirac pointsi^ A general quartic term 
allowed by the lattice symmetry, on the other hand, has 
a much smaller symmetry, as already mentioned. Nev- 
ertheless, we will argue that near the metal-insulator 
quantum critical points, all, or nearly all of the larger 
symmetry possessed by the non-interacting part of the 
Lagrangian gets restored. This conclusion is supported 
by the, admittedly uncontrolled, but nevertheless quite 
informative one- loop calculation. First, we find three dis- 
tinct critical points in the theory, all of which have not 
only the rotational, but the full Lorentz-symmetric form. 
This immediately implies that the dynamical critical ex- 
ponent is always z = 1. This is quite remarkable in light 
of the fact that the microscopic theory is not even rota- 
tionally invariant, and that the critical points in question 
are purely short-ranged.^^ The fact that z — 1 has im- 
portant implications for several key physical observables 
near the critical point, as we discuss shortly. Further- 
more, we find that two out of three critical points in the 
theory exhibit a full chiral symmetry as well. We iden- 
tify the three fixed points in the theory as corresponding 
to three possible order parameters, or "masses" that de- 
velop in the insulating phase at strong coupling. 

1) (^'7355'), which preserves chiral, but breaks time- 
reversal symmetry. Microscopically, this order parameter 
may be understood as a specific pattern of circulating 
currents, as discussed in the pasti^ 

2) (5'^'), which preserves the time-reversal symmetry, 
and the single chiral generator 735 , which will be shown 
to correspond to translational invariance. This order pa- 
rameter describes a finite staggered density, i. e. the 
difference between the average densities on the two sub- 
lattices of the honeycomb lattice,= 

3) (5* (73 cos a -f 75 sin a)\E') , which preserves the time- 
reversal, but breaks translational invariance (735). This 
order parameter can be understood as the specific 
"Kekule" modulation of the nearest-neighbor hopping 
integrals. 

In one- loop calculation all three critical points have the 
same correlation length exponent — 1, which we believe 
is an artifact of the quadratic approximation. The result 
that the dynamical critical exponent z = 1 is, on the 
other hand, possibly exact. If we denote the relevant 
interaction parameter with V, the Fermi velocity near 
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FIG. 2: The behavior of the the Fermi velocity(t;F), strength 
of the quasiparticle pole (Z), and the gap (m) near the metal- 
insulator transition. 

the transition scales as 

VF ~ [V, - vr^'-'^ (3) 

so the above value of z would simply imply that it stays 
regular at the transition. This appears to be in agreement 
with the picture of the transition as the opening of the 
relativistic " mass" in the spectrum. The mass-gap in the 
insulating phase scales as usual^^ as 

m^iV-VcYr (4) 

The transition on the metallic side is manifested as van- 
ishing of the residue of the quasiparticle pole^^ 

Z^{Vc~Vr^^. (5) 

where we assumed z ~ I. (A more general power-law is 
discussed in the text.) At one- loop the fermion anoma- 
lous dimension rj^j, vanishes, but in general it is a positive, 
small, and critical-point-dependent number. The overall 
picture of the metal-insulator transition that emerges is 
presented in Fig. 2. 

For graphene's Pz-orbitals well localized on carbon 
sites, a further significant simplification takes place. All 
the terms without the equal number of creation and an- 
nihilation operators for each of the two sublattices must 
vanish. Assuming again the emergent Lorentz symme- 
try at low energies this allows one to finally write the 
simplest internally consistent interacting theory as 

L = Lo +5D2(*735*)^ +5ci(*^')^- (6) 

This Lagrangian provides the minimal low-energy de- 
scription of interacting spinless electrons on honeycomb 
lattice. It has two critical points, corresponding to tran- 
sitions into insulators 1) and 2) in the above, both cor- 
responding to the Gross-Neveu criticality in 2-1-1 dimen- 
sions. We discuss the internal consistency and the suf- 
ficiency of this Lagrangian and some of the peculiarities 
of the ensuing phase diagram. 



The rest of the paper is organized as follows. We dis- 
cuss the point symmetry, translational symmetry, and 
the time-reversal symmetry of the interacting Lagrangian 
as dictated by the microscopic Hamiltonian for the sys- 
tem in the next section. In section III it is shown how 
further enlargements of the symmetry would reduce the 
number of coupling constants. We introduce the notion 
of "maximally symmetric" theory, which shares the full 
Lorentz and chiral symmetry with the quadratic term in 
the Lagrangian. The general formalism of Fierz transfor- 
mations is developed and applied to the cases of interest 
in section IV. The change of the coupling constants with 
the ultraviolet cutoff in the theory is studied in section 
V. The atomic limit of the general interacting theory is 
described in section VI, and the critical exponents are 
discussed in section VII. In section VIII we discuss the 
scaling of the electron propagator and the power-laws for 
various quantities of interest. The discussion of the long- 
range Coulomb interaction and the critical behavior of 
the specific heat and the optical conductivity are given 
in section IX, and the summary in section X. Finally, 
in Appendixes we present some of the requisite techni- 
cal details: the Fierz transformation, the spectral form 
of the asymmetric matrix needed in section IV, and an 
alternative implementation of the renormalization group 
in presence of linear constraints. 

II. SYMMETRIES AND SHORT-RANGE 
INTERACTIONS 

A. Hamiltonian and the Lagrangian 

As the simplest microscopic model that contains 
the relevant physics we may consider the tight-binding 
Hamiltonian on the graphene's honeycomb lattice, de- 
fined as 

Ho = t^u\A)viA + b',) + H.c., (7) 

A,i 

where u and v are the electron annihilation operators 
at two triangular sublattices of the honeycomb lattice. 
Here, A denotes sites of the sublattice generated by 
linear combinations of basis vectors di = (a/S, — l)a, 
02 — (0, l)a, whereas B = A + b are the sites on the 
second sublattice, with 6 being bi — (l/\/3, l)a/2, 62 = 
(l/V3,-l)a/2, or bs = (-l/%/3,0)a, and a is the lattice 
spacing. 

Within the framework of the tight-binding model 
the energy spectrum is doubly degenerate E{k) = 
±t|^jexp[/c • bi\\, and becomes linear and isotropic in 
the vicinity of six Dirac points, at the edge of the Bril- 
louin zone, among which only two, hereafter chosen to be 
at ±K with K = (1, l/-\/3)(27r/a'\/3), are inequivalent. 
Retaining only the Fourier components in the vicinity 
of these two inequivalent points, the quantum mechan- 
ical action corresponding to Hq at low energies can be 
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written in the form 5 = /q drdxLo, with the free La- 
grangian density Lq defined as in Eq. (2), with r as the 
imaginary time and T is the temperature. Matrices 7^ 
satisfy the Chfford algebra {7^, 7^} — 2Sf^^, /i, = 0, 1, 2, 



Here, the reference frame is conveniently rotated so that 
= q- K/K, qy = iKxq)x K / , Lu,, = (2n + 1)t:T 
are the fermionic Matsubara frequencies, A 1/a is a 
high-energy cutoff, and we set h = ks = vp = ^, where 
vp = ta a/3/2 is the Fermi velocity. Choosing 

implies 

The two remaining anticommuting matrices may then be 
taken as 

7,^ (:/„'), 7.^(4--')^ (10) 

This defines the "graphene representation" of 7- 
matriceS)^ a are the standard Pauli matrices. 

Considering a more general model with further hop- 
pings or weak anisotropics^ can be seen not to destroy 
the Dirac points, but only to shift them in energy. As 
this can always be compensated by a shift of the chemi- 
cal potential, the Lagrangian (2) provides the low-energy 
description of the general free electronic Hamiltonian on 
a honeycomb lattice, with the chemical potential tuned 
to the Dirac point. 

Note that the free Lagrangian besides the Lorentz sym- 
metrt also possesses another chiral, pseudospin or "val- 
ley" , global SUc{2) symmetry, generated by {73, 75, 735}- 
Both the Lorentz and the chiral symmetry of the free 
Lagrangian are emerging only at low energies, however, 
and the term quadratic in derivatives in Lq, for example, 
would spoil it. As will be shown shortly, both symme- 
tries are also violated by the leading irrelevant quartic 
couplings introduced by the interactions. 

Let us now consider the electron-electron interactions. 
The Hamiltonian of a general four-fermion interaction 
has the form 

H^nt^ ^ {af3\V\^5)r^\0^rsr^, (11) 



^ = ^'l'7g. The summation over repeated space-time in- 
dices is assumed hereafter. The fermionic field '^'{x, t) is 
defined as 



,v\K + q, ujn), u\-K + q, u,,), -yt {-K + q, w„)). (8) 

I 

where r — u ov v are fermionic annihilation operators, 
and the matrix element corresponding to the interaction 
potential V{r) is given by 

{al3\V \ -fS) = J dxdyipl{x)ip)i{y)V{x - y)ip^{x)fsiy)- 

(12) 

Here we can take fa{x) to be a localized px-orbital on 
the site a, so that it belongs to either one of the two 
sublattices of the honeycomb lattice. In general, there 
is no restriction on the overlap of the wave-functions, 
and all the matrix elements (a/3 | V \ ^5) are in principle 
finite. Their relative sizes, however, may be rather differ- 
ent, and we discuss important simplifications that follow 
in the limit of well localized orbitals in sec. VI. In the 
following we will consider general "short-ranged" inter- 
actions, which are defined by the interaction V{x) with a 
regular Fourier component V{k — 0). Without a loss in 
generality one may then take the interacting Lagrangian 
for spinless fermions corresponding to Hint as in Eq. (1) 
where Mi and M2 are some constant 4x4 Hermitian 
matrices. The interacting Lagrangian contains therefore 
at most 16-1- (16 x 15/2) — 136 independent real coupling 
constants. However, the number of couplings in Lint is 
severely reduced by the lattice symmetries, as we discuss 
next. 



B. Reflection symmetries 

Two obvious discreet symmetries of the honeycomb 
lattice that have not been broken by our choice of the 
Dirac points are the reflection symmetries through the 
lines A and B in Fig. 1. Let us consider the former sym- 
metry first. It exchanges the two sublattices, but not the 
two Dirac points. The low-energy Lagrangian thus has 
to be invariant under the exchange of the spinor com- 
ponents belonging to different sublattices, u{k)^v{k). 
Consequently, the symmetry operator acting on the four- 
component Dirac spinor defined in Eq. ([H]) has the form 

S^l®a^^-12. (13) 

Since under this reflection q^ — > qx and qy — > — is 
evidently invariant under S. The invariance of Lint under 
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this reflection symmetry requires both matrices Mi and 
M2 to either commute or anticommute with the operator 
S: 

[S,Mi] = [S,M2]=0, (14) 

or 

{S,Mi} = {S,M2} = 0. (15) 

Similarly, the reflection symmetry through the line B 
exchanges the two Dirac points, while not exchanging the 
sublattice labels. It corresponds therefore to 

^ = ^^^^^ = (/, o)- (16) 

Recalling that under this transformation —Qx 

and Qy Qy, it is evident that the free Lagrangian in 
graphene representation is indeed invariant under T as 
well. Demanding the interacting Lagrangian Lint to be 
invariant under the action of the operator T on the Dirac 
spinor, implies that both matrices Mi and M2 either 
commute or anticommute with T as well. In other words, 
both matrices Mi and M2 have to be either even or odd 



with respect to T: 

[T,Mi] = [T,M2]=0, (17) 

or 

{r,Mi} = {r,M2} = o. (18) 

Together with the combination of the two reflections S 
and T, and the identity operation, the two symmetry op- 
erations form the dihedral group (or Klein's vierergrouppe, 
in older literature) D2: D2 = {1,S,T,ST} = Z2 x Z2, 
the symmetry group of a rectangle. Note that the trans- 
formation ST is just the space inversion, and that the 
rotation by tt/2 does not belong to D2. 

One may now classify all the four-dimensional ma- 
trices into four categories, according to their trans- 
formation under the two reflection symmetries S and 
T, respectively: even-even, A = 72, i7o73, ^7175}) 
even-odd, B = {17071^ 735, no75, ^7173}, odd-even, 
C = {7o,«7o72,73>n273}> and odd-odd, D = 
{71, Z7i72, 75, «7275}- The interacting Lagrangian sym- 
metric under the D2 is thus restricted to be of the fol- 
lowing form 



Lint = aij{^^Ai^){^^Aj^) + bij{^^Bi^){^^Bj^) + Cij{^^Ci^){^^Cj^) + rfi,(*t£).^-)(^-t£),^'), (19) 

I 



where o,j, o = a,b,c,d, i,j = 1,...,4, are real and sym- 
metric. The maximal number of independent real pa- 
rameters specifying the allowed couplings is thus already 
reduced to forty, since each Oij has ten independent com- 
ponents. 

C. Translational invariance 

The generator 735 = az ® h of the chiral symmetry 
plays a special role. It is in fact the generator of trans- 
lations. To see this recall that under a translation by R 
the electron fields transform as 

r(fc,w) ^ e'^-^r(fc,w) (20) 

where r = u,v. The Dirac field under the same transfor- 
mation thus changes as 

^{q,u) ^ e'^^-^^''''e*^-^^{q,ij), (21) 
or, in real space, 

*(f,r)^e*(^-^)^^^*(£+i?,r). (22) 

Translational invariance requires therefore that Lint is a 
scalar under the transformations generated by 735 , which 



we will denote asUc{l). It is easy to see that this is pre- 
cisely the same as demanding the conservation of momen- 
tum in the interaction terms. The reader is also invited 
to convince herself that the terms with the higher-order 
derivatives in Lq would also be invariant under the ?7c(l)- 
First, we observe that there are eight linearly indepen- 
dent bilinears that are scalars under the Uc{l): 

Xpi = (23) 

where F = A,B,C,D and i = 1,2. The remaining eight 
bilinears can be grouped into four vectors under the same 
C/c(l): 

a = (*U34',*^53*), (24) 
/3 = (*tB4^,^U4^), (25) 
7 = (*'fC3*,*tD3^»), (26) 

5= {<i/^C4^,^^D4^), (27) 

The invariance under f/c(l) implies therefore that the 
interacting Lagrangian has the following form 

Lint = ^ gPiXpi + ^ gFXFlXF2 (28) 
Fi F 
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The number of possible independent couplings is down 
to eighteen. 

D. Time-reversal 

The set of the allowed couplings is further reduced by 
the time-reversal symmetry. The microscopic interact- 
ing Hamiltonian (jlip is invariant under the time-reversal, 
and therefore the corresponding low-energy Lagrangian 
has to possess the same invariance. The time-reversal 
symmetry requires that ItHIt~^ = H, where It is the 
antiunitary operator representing the time-reversal sym- 
metry, and thus has the form It = UK, with U rep- 
resenting the unitary part of It and K is the complex 
conjugation. To find the form of It let us consider first 
the massive Dirac Hamiltonian 

H = ijojiPt + mi7o, (29) 

with the mass mi describing the imbalance in the chemi- 
cal potential on the two sublatticesi^ Recalling that mo- 
mentum changes sign under the time- reversal, ItPiIt~^ = 
—Pi, in the graphene representation the invariance of the 
above Hamiltonian under the same transformation im- 
plies 

{[/, ijoli} = [U, Z7072] = [U, 7o] = 0, (30) 

and hence C/ = ie*'^7i(cos6'73-|-sin6'75). Within the sim- 
plest framework of the tight-binding model with uniform 
hopping the time-reversal operator is not uniquely de- 
termined. We thus consider a generalized tight-binding 
model with anisotropic hopping defined as 

Ham.o = J2^i+di^^y{A)v{A + h) + H.C., (31) 

A,i 

where 

= ^-/^{A)e'^'^^e'^-^ + c.c. (32) 

represents a non-uniform hopping, and G = 2K M- On 
a lattice, this particular set of hoppings generates the 
so-called Kekule texture. Near the two Dirac points the 
Hamiltonian Haniso reads 

Haniso = holiPi + ™2«7075 + ^,317073, (33) 

where m2 = /m(A(r)) and = i?e(A(r)). The two 
masses TO2 and therefore provide the low-energy rep- 
resentation of a completely real microscopic Hamiltonian, 
so that we postulate that Haniso is also time-reversal 
symmetric. In graphene representation this requires the 
unitary part of the time-reversal operator to obey the 
following algebra: 

[C/,i7o73] = {(7,17075} = 0. (34) 



The matrix T = 17175 satisfies conditions (p4|) and 
thus the unitary part of the operator It acting on the 
spinless Dirac field ([8]) is^^ 

U = T = Z7175 =(^l ) , (35) 

with I2 as the 2x2 unity matrix. The unitary part of the 
time-reversal operator thus simply exchanges the compo- 
nents of the Dirac spinor "^{x, r) with different valley in- 
dices, as expected. It also happens to be the same as one 
of the two matrices representing the reflection operators 
from D2. 

Another way of arriving at the same form for the time- 
reversal operator is to postulate that an arbitrary chi- 
ral transformation of the Dirac Hamiltonian in Eq. (29) 
yields a time-reversal invariant Hamiltonian. Alterna- 
tively, our derivation may be understood as a demon- 
stration of commutativity of the chiral and time-reversal 
transformations. 

Since we have already used the invariance under T to 
restrict the interacting Lagrangian, time-reversal invari- 
ance will be observed if the remaining terms are even 
under complex conjugation. All the terms Xp^ and 
p- p are thus automatically invariant under time- reversal, 
but among the remaining six mixed terms, the terms 
XciXc2, X£)iXd2, and 7 • 6 are odd. Time-reversal 
invariance implies therefore that 

9c = 9D= 9'tS = 0, (36) 
which leaves then at most fifteen independent couplings. 



III. ENLARGEMENT OF SYMMETRY 

We found that the exact symmetries of the microscopic 
Hamiltonian, D2, translational, and the time-reversal, 
leave at most fifteen independent short-range couplings. 
Anticipating some of the results, it is interesting to de- 
duce the further reductions of the number of couplings if 
one by hand imposes larger symmetries onto the interac- 
tion Lagrangian Lint ■ 



A. Rotational invariance 

Since the rotation by 7r/2 is not a member of the Z?2, 
the matrices 71 and 72 appear asymmetrically in Lint- 
If we demand that they appear symmetrically. Lint be- 
comes fully rotationally invariant. This is achieved if 

gA = 9b = gap = gA2 - goi = gsi - gc2 = gp- gs = 0. 

.(37) 

Let us call the interacting Lagrangian with the rotational 
invariance imposed this way Lint.rot- It would be de- 
scribed by at most nine couplings. 
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B. Lorentz invariance 



A. General problem 



Imposing further the Lorentz invariance would require 
that on top of the above restrictions one also has that 



gAi + gsi = gA2 + 932 = g-r + gf3 = 0. 



(38) 



With the restrictions in the previous two equations the 
Lagrangian has only six couplings constants, and may 
be cast in a manifestly Lorentz invariant form, worth 
displaying: 

^int,lor — gAl 

(i*75*)^] +57[(*7m73*)^ + (%^75*)^]- 

C. Chiral invariance 

Finally, the maximally invariant interacting La- 
grangian would be with the full, i. e. both the Lorentz 
and the chiral, symmetry of the non-interacting part. 
This is achieved by setting in the last equation 



gci - ga = gB2 - = 0. 



(40) 



The interacting Lagrangian can in this case be written 
as 



Lmt,max — g AlS^^ + gD2S^ + gciV + 552^^^: 

where the participating bilinears in Dirac fields. 



*7m*. 



S = *735*, 



V = (\I'*,i\E'73\I',i\E'75\I'), 



= (*7m735*, *7m73*, *7m75*), 



(41) 

(42) 
(43) 
(44) 
(45) 



are the scalar (vector), scalar (scalar), vector (scalar), 
and vector (vector) under the chiral (Lorentz) transfor- 
mation. The last form makes the Lorentz and the chiral 
symmetry of the Lagrangian LQ+Lint,max manifest. Such 
a maximally symmetric Lagrangian contains therefore at 
most only four coupling constants. 



The application of Fierz identity to the set of quar- 
tic terms allowed by the assumed symmetry in principle 
leads to the set of linear constraints of the form 



FX = 0, 



(46) 



where is a real typically asymmetric matrix, and X is 
a column; the elements of which are the quartic terms al- 
lowed by the symmetry. Of course, only the quartic terms 
which share the same symmetry may be related by Fierz 
transformations. For example, in the maximally sym- 
metric case = (5'^, 5^, T?^, l7f ). When the number 
of couplings is small it is easy to discern the linearly in- 
dependent combinations of the original terms, but when 
it is not, as the case is for D2 x ?7c(l) x It microscopic 
symmetry of Lint, one needs a more general method of 
doing so. 

In Appendix B we show that an asymmetric matrix 
such as the Fierz matrix F can be written in the diadic 
form 



29 



(47) 



where {fJ-i} is the real spectrum of the related symmetric 
matrix Sp — F^ F. In the eigenbasis of Sp we can write, 
in Dirac notation. 



(48) 



so that the above linear equations can be written as 



0. 



(49) 



Since the vectors also form a basis, it must be that 

either: 

a) for fii ^ 0, {fJ.i\X) = , or 

b) = 0, so that {^J.^\X) ^ 0. 

The first set provides us then with the linearly indepen- 
dent constraints, and the second with the set of remaining 
linearly independent quartic terms. Since the matrices F 
and Sp obviously have the same kernels, the number of 
independent coupling constants allowed by the symmetry 
is simply the dimension of the kernel of the appropriate 
Fierz matrix. 



IV. FIERZ TRANSFORMATIONS 

The number of independent couplings is further re- 
duced by the existence of algebraic identities between 
seemingly different quartic terms. The derivation of the 
so-called Fierz transformation, which allows one to write 
a given local quartic term in terms of other quartic terms 
is provided in Appendix A. A systematic application of 
this transformation allows one to reduce the number of 
independent couplings for a given symmetry of the inter- 
acting Lagrangian. 



B. Maximally symmetric case 

Let us consider the simplest example of the quartic 
term with the full Lorentz and chiral symmetry, Lint,max 
first. Defining the vector X as in the above leads to the 
Fierz matrix 



F 



/ 3 3 3 -1 
5 111 
3 3 3 -1 

V9 -3 -3 5 



(50) 
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with the two-dimensional kernel with the zero-eigenstates 



D. Lower symmetries 



1 



(/ill = -^(0,-1,1,0), 



(M2| = ^(-1,1,1,3), 



(51) 



(52) 



and Hi = fi2 = 0. The remaining two eigenvalues are 
/Z3 = 64, and fi^ — 144. The general method explained 
above implies that the general maximally symmetric in- 
teracting Lagrangian can be written as 

L^nt,macc = Xl{V^ ^ Sf,) + X2{- + S^ + V^ +3V^) , (53) 

with the "physical" couplings 

gci - 9A1 



A. = 



Ai 



-9D2 + gAl + gCl + 3ffB2 

12 



(54) 



(55) 



The two remaining linearly independent combinations 
vanish due to Fierz identity. So the maximally symmetric 
interacting theory is specified by only two quartic cou- 
pling constants, which may be chosen to be any linearly 
independent combinations of the above Ai and A2. 



C. Lorentz-symmetric case 



The symmetry ladder may now be climbed back to 
the £>2 X C^c(l) X It minimally symmetric Lagrangian. 
First, reducing Lorentz to rotational symmetry increases 
the number of independent couplings to four. Remov- 
ing the rotational symmetry finally increases the number 
of couplings to six. We may note in passing that there 
are two independent Fierz identities between the mixed 
terms that obey the time- reversal symmetry in Eq. (28): 

?,XaiXa2 - XbiXb2 + dx0^Q, (61) 



XaiXa2 + XbiXb2 + a X /3 = 0, (62) 

so that the three mixed terms in fact contribute a single 
independent coupling. 



V. RENORMALIZATION GROUP 

Having determined the independent coupling constants 
for each symmetry, we now proceed to study their 
changes with the decrease of the upper cutoff A. We will 
be particularly interested in fixed points of such renor- 
malization group transformations, as they will provide 
the information on the quantum metal-insulator transi- 
tions that can be induced by increase in interactions. 



We may then proceed to find the independent cou- 
plings for the next case in order in complexity, Lint.ion 
with the chiral symmetry broken down to t/c(l)- If we 
define the six-dimensional vector 



2 ^ 



X ' = (S^ sl,Vi^ vi + vi, s^, sl^ + s. 

the Fierz matrix is found to be 



F 



with the three-dimensional kernel spanned by 
1 



/I 


1 


5 


-1 


1 


-1 \ 


3 


3 


-1 


5 


-7/3 


-1/3 


3 


3 


3 


3 


-1 


-1 


5 


1 


1 


1 


1 


1 


3 


-1 


3 


-3 


3 


1 


V3 


-1 


-3 





1 


2 / 



5%/2 



(-2,3,1,0,0,6), 



1 



5V2 
1 



2V55 



(-1,4,-2,0,5,2), 
(-3,-7,3,10,7,2). 



(56) 



(57) 



(58) 



(59) 



(60) 



Lint.ior can now be written in terms of only three linearly 
independent quartic terms, in complete analogy with the 
maximally symmetric case. 



A. Maximally symmetric theory 

Let us again begin with the maximally symmetric La- 
grangian, L — Lo+Lint.max- There are only two coupling 
constants to consider in this case, and we choose them to 
be gD2 and qai, which correspond to S'^ and quartic 
terms, respectively. If any of the other two terms would 
become generated by the renormalization transformation 
we would use the Fierz identity to rewrite it in terms of 
and S^. Alternatively, one may wish to renormal- 
ize the theory as written in terms of physical couplings 
in Eq. (53). This procedure, completely equivalent to 
what is pursued here, is described in Appendix C. As we 
integrate the fermionic modes lying in the 2+1 dimen- 
sional momentum shell'''' from A/b to A, with 6 > 1, to 
quadratic order in coupling constants we find 



dg 



D2 



d\nb 



-9D2- gl,2 + '^9Al+^9D29Al, (63) 



= -9A1 + 9ai + 9029 AI- (64) 

We rescaled the couplings here as 2gA/TT'^ — > g. To this 
order no other types of quartic terms actually get gener- 
ated, and the Fierz transformation turns out not to be 
necessary. The limit of the above equations that survives 
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FIG. 3: Schematic flow diagram of the two coupling constants 
in the maximally symmetric theory. The flxed point A de- 
scribes the continuous transition into a time-reversal symme- 
try broken insulator, and C the dynamical generation of the 
chiral symmetry breaking mass. The line gD2 = describes 
the Thirring model, and the dashed line the Nambu-Jona- 
Lasinio model in 2-1-1 dimensions. (See the text.) 

the extension to a large number of Dirac fields also agrees 
with the previous calculationi^ 

The above flows, besides the Gaussian, exhibit three 
fixed points at finite couplings (Fig. 3) . The first critical 
point (A) is at 5/32 — .9^11 — 0, and the second critical 
point (C) is at gD2 = (^5 - l)/2, gAi = (3 - V5)/2. 
There is also a bicritical fixed point (B) that separates 
the domains of attraction of the two critical points, at 
gD2 = -{V5 + l)/2, gAi = {V5 + 3)/2. 

The physical interpretation of the critical point A is 
obvious. Since we can tune through it by keeping gAi = 
and increasing g]j2 over a certain negative value, it should 
describe the transition into the insulator with the gap 
that breaks the time-reversal symmetry, described by 

(^ss^-) ^ 0. (65) 

This state is obviously favored at a large and negative 
gD2- Note that since 735 commutes with 7^, the line 
gAi = is invariant under RG. In fact, the perturba- 
tive /3-function along this line has to be identical as the 
one in the Gross-Neveu model. We can therefore simply 
use the already existing higher-order estimate a'^^i'^^ and 
the numerical results'^'''^^ to find the critical exponents 
describing this particular metal-insulator transition. We 
return to this fixed point shortly. 

The physical interpretation of the critical point C is 
less obvious, but we can think of it as follows. First, 
note that the line gD2 = 0, gAi > 0, which describes the 
Thirring model,- ^ belongs to the domain of attraction of 
C. Also, the Fierz transformations in Eq. (50) imply 

gD2S^ + gAiSl = {gD2 - '2gAi)S^ -~ gAiV^- (66) 



The line gD2 — '^gAi = for gAi > , which we name the 
Nambu-Jona-Lasinio (NJL) liner^ also falls into the do- 
main of attraction of the critical point C. Along this line, 
however, there should be a transition into an insulating 
state with 

(n-V)^ 0, (67) 

where n is a unit vector. Such a state is clearly favored at 
a large and positive gAi along the NJL line, and breaks 
the chiral SUc{2) symmetry down to U{1). We therefore 
identify C as the metal-insulator critical point governing 
the chiral-symmetry breaking transition in both Thirring 
and NJL models with a single Dirac field. 

The picture suggested by the above one-loop calcula- 
tion in the maximally symmetric theory appears quite 
natural. There are two possible insulating phases, each 
breaking either chiral or time-reversal symmetry, which 
correspond to possible "masses" for the Dirac fermions. 
Both metal-insulator transitions are continuous, and are 
described by different critical points. 

Of course, the true low-energy theory on the honey- 
comb lattice is much less symmetric than the one studied 
in this section. Nevertheless, we will argue that the two 
identified critical points may in fact be stable at least 
with respect to weak manifest breaking of the Lorentz 
and chiral symmetries. 



B. Broken Lorentz symmetry 

The explicit breaking of Lorentz symmetry down to 
the rotational symmetry can be easily implemented by 
adding to Lint,max a small symmetry breaking term 

<5(*7o«')', (68) 

which by virtue of being only rotational symmetric is 
guaranteed to be linearly independent of S'^ and S^. A 
weak perturbation 5 to the lowest order in couplings 
and gAi then fiows according to 

■^^=5{-l-gD2+gAi) + 0{5^). (69) 

We thus find the critical point C to be stable with re- 
spect to weak Lorentz symmetry breaking to one loop, 
the bicritical point B to be unstable, and A marginal. 
We suspect that this result, although clearly an outcome 
here of an uncontrolled approximation, may be indica- 
tive of the true state of affairs. Hereafter we will as- 
sume that the critical points A and C are stable with 
respect to weak breaking of the Lorentz symmetry in the 
Lagrangian. It may also be worth mentioning that the 
complete one- loop /3-functions for 5, gAi and gD2i which 
we have computed but have not shown, do not lead to 
any new critical points at 5 7^ 0. 
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C. Broken chiral symmetry 

The simplest quartic term with the full Lorentz sym- 
metry, and only Uc{l) subgroup of the full chiral symme- 
try may be written as 



Lvnt,lor = gD2S^ + gciV^ + gaiV^ 



V4 



(70) 



The Fierz transformation matrix given above implies that 
these three quartic terms are indeed linearly independent. 
When gc\ — gen the Lagrangian Lint.ior acquires the 
full chiral SUc{2) symmetry, and may be rewritten as 

Using the Fierz transformation, and after a conve- 
nient rescaling of the couplings as gA/Sn^ — > g, to the 
quadratic order one finds 



dg 



D2 



d\nb 



-gD2 - 6gD2 



(71) 



45a + 6gD25ci + 12gn25a - ^gciga. 



dgci 

d\nb 



-gci - &gci - ^gl + 6gu25ci - 4:gciga, (72) 



dgg 
d\nb 



la + ^gD2ga " lOgagCl- 



(73) 



The two chirally symmetric critical points from the 
previous section now appear at gD2 — —1/6, gci = ga — 
(A), and gD2 = (3^5 - 7)/12, gci = ffa = (VS - 
3)/12 (C), and both remain critical, even in absence of 
chiral symmetry in the Lagrangian. There is, however, 
an additional critical point (E) at gjj2 = (\/5 — 2)/6, 
gci = — 2(7q. — (-\/5 — 3)/6. Note also that the plane 
(7c( = is invariant under the renormalization group, but 
whereas the fixed point A in that plane is critical, the 
fixed point (D) at — 0, gci = —1/6 is bicritical. One 
also finds that the line gD2 = 0, 1701 < 0, ga-infinitesimal 
and positive intersects the critical surface which contains 
the point E, whereas for i^Q-infinitesimal and negative the 
critical behavior is governed by C. For weak ga therefore 
there is a crossover from the fixed point at D toward 
either C or E, depending on the sign. Interestingly, for 
negative ga chiral symmetry becomes fully restored at 
the transition, at least within our one-loop calculation.^^ 



VI. ATOMIC LIMIT 

Motivated by the one-loop results, we will assume here- 
after that the Lorentz symmetry becomes restored at long 
distances in the domain of interest, and that we need only 
consider Lint.ior with the three couplings from the last 
section. The situation however, can then be simplified 
even further, as we discuss in this section. 

Consider the interaction Hamiltonian in Eq. (11). If 
the p^-orbitals are well localized on their corresponding 



lattice sites, we may neglect the matrix elements with 
a 7^ 7 or /3 ^ 5. Keeping only the remaining, dominant 
matrix elements then one obtains the "atomic limit" of 
the general interaction Hamiltonian 



Hi, 



Hi 



lat 



E 

a,l3 



Va^prianp 



(74) 



where Ua is the electron number operator at site a. 
The class of Hamiltonians Hiat is evidently still rather 
broad, and would for example include all lattice interact- 
ing Hamiltonians. 

Writing the lattice Hamiltonian Hiat in terms of the 
Dirac fields, however, imposes yet another restriction on 
the coupling constants. Since Hiat is written in terms of 
lattice-site particle number operators, any Dirac quartic 
term evidently must contain an equal number of (w^) 
and u (v) fields. On the other hand, the g^-term from 
above in momentum space can be written schematically 
as 



{Vi + V^) - {u\v2 + vlu2){ulvi + vlui) 



(75) 



and thus contains the terms forbidden in the atomic 
limit, "^^ such as u\v2U2Vi. The index 1 and 2 refers here 
to the two Dirac points. This implies that for any lattice 
Hamiltonian Hiat we must have 



ga = 0. 



(76) 



Note that the plane ga = is invariant under the change 
of cutoff in the above one-loop calculation. It is easy to 
see that this feature of the /9-functions for gD2, gci and 
ga is in fact true to all orders in perturbation theory. The 
matrices 735 and / in the remaining two terms in Lint.ior 
commute with the Dirac propagator, and therefore an ar- 
bitrary diagram containing g]j2 and gci terms can con- 
tribute only to the renormalized go2 and gci couplings. 
So imposing ga — at an arbitrary cutoff guarantees its 
vanishing at all others. 

It is therefore not only physically justified but also in- 
ternally consistent to consider only the two couplings go2 
and gci in the Lorentz symmetric, but chirally asymmet- 
ric low-energy theory. The one-loop result in this plane 
is depicted in Fig. 4. The transition is either into the 
time-rcversal-symmetry-broken, or into chiral-symmetry- 
broken insulator. Several features of this ffow diagram 
that should be generally valid are worth mentioning. 

1) There should be two critical points, both unstable in 
a single direction. Bicriticality of A, for example, would 
imply that the transition for negative gD2 at a weak pos- 
itive gci is first order. This, however, seems unlikely on 
physical grounds, and, also, it is not found in the ex- 
plicit large-N generalization of the theory, when the two 
/3-functions are known to decouplei^ 

2) The two critical points have identical critical be- 
havior. This is because the term ^'735 5*, in graphene 
representation, under the transformation 

* ^W2 + cTz) (8) tj^ + {I2 - <Jz) ® h]^, (77) 
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FIG. 4: The flow diagram in the gD2 — gci plane. The possible 
non-perturbative fixed point at the gD2 = gci line, which gov- 
erns the transition in the two-component Gross-Neveu theory 
is not shown. 



First, since all the identified critical points exhibit 
Lorentz symmetry, 

z = l. (81) 

We also find that the exponent v is unity at all criti- 
cal points as well, but this is clearly an artifact of the 
one-loop calculation. In general, v is expected to be dif- 
ferent at different critical points. The same goes for 77^, 
which vanishes in one-loop calculation, but will be finite 
in general. 

In the atomic limit, when under the assumed Lorentz 
invariance we need only two coupling constants, the val- 
ues of the critical exponents are better known. First, 
in a perturbative calculation in powers of coupling con- 
stants, the exponents at critical points A and D will be 
identical. We may thus expect that in a lattice theory 
with short-range repulsion the transition is either into 
the time-reversal symmetry or chiral symmetry broken 
insulator, in cither case with^ 

1^ = 0.74- 0.93 (82) 



^rt ^ l*t[i(/2 + a,) ®a, + {h - (Jz) ® I2], (78) 

goes into ^5", and vice versa, while Lq remains invariant. 
This also means that the two /3-functions are symmetric 
under the exchange gD2 gci- Both critical points are 
thus in the universality class of the Gross-Neveu model. 
3) At the line §02 — gci the single /^-function becomes 



dg 



D2 



d\nb 



(79) 



i. e. gD2 flows according solely to its canonical dimen- 
sion. This is because 



(^-735*)" + (**)' = 2{^\az^+Y + 2{^laz^-Y, (80) 

where = (^'^,5'^). Since all 7^ are block-diagonal, 

~\-K and —K components at this line decouple. The par- 
tition function factorizes into a product of two Gross- 
Neveu partition functions, each containing a single two- 
component Dirac fermion. Along this line the system is 
believed to have the metal-insulator transition, possibly 
continuous, but the /3-function vanishes at least to the 
order gjj2'^ 



VII. CRITICAL EXPONENTS 

Each identified metal-insulator transition is character- 
ized by a set of critical exponents. We will here focus 
on the three already mentioned in the introduction: the 
correlation length exponent v, the dynamical exponent 
z, and the Dirac fermion anomalous dimension 77^. The 
other exponents can then be obtained from the usual 
scaling relations)^ 



r]^ = 0.071 -0.105. 



(83) 



Note that since there is only a single Dirac field involved 
the numerical values of the exponents v and rj^j, differ 
significantly from the large-N values of unity and zero, 
respectively. This raises hope that this non-trivial critical 
behavior may be observable in numerical simulations of 
lattice models. 

Although not of immediate relevance to graphene, it 
would still be of interest to determine the critical expo- 
nents at the chirally symmetric critical point C, which we 
proposed to control the critical behavior of the Thirring 
and the NJL model. We, however, are not aware of any 
analytical nor numerical study of the NJL model that 
goes beyond the leading order in large-N calculation4£ 



VIII. FERMI VELOCITY AND RESIDUE OF 
QUASIPARTICLE POLE 

The critical exponents, as usual, govern for example 
the critical behavior of the gap on the insulating side 
of the transition, as mentioned in the introduction. In 
the present case, however, there are massless fermionic 
excitations on the metallic side, and one may wonder if 
and how the approach to the critical point is refiected 
onto these. Let us therefore generalize slightly and pro- 
vide the support for the results already announced in the 
introduction. 

First, the usual scahngS^ imphes that at the cutoff A/6 
electron's two-point correlation function near the critical 
point and at zero temperature satisfies 



G = h''F{hk,b'uj,tb^'''), 



(84) 
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where oj is the Matsubara frequency and t ~ (14 — > 
is the transition's tuning parameter. Setting tb^/^ — 1 
we thus find the usual scahng law 



G = t-'"'F{t-''k,t-'''ijj) 



(85) 



where F{xi,X2) = F(xi, X2, 1) is a universal scaling func- 
tion. From here we can extract the scaling of the Fermi 
velocity and quasiparticlc residue as follows. First, if 
upon the analytical continuation to real frequencies G 
has a pole at 

uj = VF{t)k (86) 

the scaling relation immediately dictates that 

VF{t) = e^'-^^F. (87) 

Let us next set w = 0, take A: > 0, and let t — > 0. In 
this limit 



G 



1 



(88) 



and therefore F{xi — > oo,0) ~ l/x\~^''. In order to 
cancel the t-dependence of the prefactor in Eq. (85) in 
this limit it must be that 



(89) 



where rjk is the (momentum) anomalous dimension. 
Analogously, assuming that for k = G ^ 
one finds that also 



z{l - Vu)- 



(90) 



On the other hand, in the opposite limit t > and 
w — *■ 0, in the metallic phase at low energies we have 
fermionic quasiparticles. This implies that, for example, 
F{0,X2 ^ 0) - 1/X2, i. e. 



G 



Z 



(91) 



with Z ^ ^^'^ . Combining with the previous relation, 
the quasiparticlc pole's residue behaves as 



(92) 



For z = 1, the two anomalous dimensions are the same, 
= ^Ilo = ^7*, where ry^ is the Dirac fermion's anomalous 
dimension, and the scaling announced in the introduction 
follows. The special form of this relation for large number 
of Dirac components when besides z = 1, it is also v — 1, 
was previously proposed.— The quasiparticlc residue van- 
ishes upon the approach to the metal-insulator transition, 
as proposed long-ago by Brinkman and Rice, but here in 
a decidedly non-mean-field fashion. 



IX. DISCUSSION 

There are at least two obvious generalizations impor- 
tant for real graphene: the addition of spin, and the 
inclusion of the long-range tail of Coulomb repulsion. 
Adding spin would simply double the number of cou- 
plings for each symmetry, since each independent quartic 
term would then require a separate coupling in the sin- 
glet and in the triplet channels. The minimal internally 
consistent low-energy theory would then be the general- 
ization of the Lorentz invariant Lagrangian in Eq. (70), 
with go, = 0: 



J spin 

intJor 



(93) 



where the sum goes over M = 1,735, and i — 0,x,y,z, 
with (To = and gn^x = gM,y = gM.z- The La- 
grangian with 5^35, i — would represent the extended 
Hubbard model with on-site and nearest neighbor repul- 
sion, considered before in the limit of large number of 
Dirac fermions in ref. 2. The interplay between the vari- 
ous instabilities in the theory equivalent to the above La- 
grangian was recently studied ir4i, where it was pointed 
out that the second-nearest-neighbor repulsion implies 
a negative coupling g-y^^^o for example. The form of the 
above minimal spinful Lagrangian facilitates a systematic 
study of the metal-insulator transition in the Hubbard 
model, which will be a subject of a separate publication. 

Few comments on the importance of long-range tail of 
Coulomb interaction are also in order. Weak Coulomb 
(~ e^/r) interaction is an (marginally) irrelevant per- 
turbation at the Gaussian fixed point, and this remains 
true at the metal-insulator critical point at large-N as 
well^i^i^. Furthermore, the entire /3-function for the 
charge coupling can be computed at large-N, and it 
does not exhibit any non-trivial zeroes.— 1^ On the other 
hand, several calculations show that by increasing the 
coupling beyond certain point and for small enough 
N the system can be tuned through a metal-insulator 
transition at which the chiral symmetry becomes spon- 
taneously broken fiii^i^i^i^ The nature of such a pu- 
tative metal-insulator transition is in our mind an open 
question at the moment. Whereas it is possible that it is 
described by new "charged" critical point^ correspond- 
ing to the non-trivial zero of the /3(e^), it also seems con- 
ceivable that the charge is always irrelevant and that the 
transition is still in the universality class of the critical 
point C, in our nomenclature. Yet another possibility is a 
discontinuous transition. More work is obviously needed 
in order to be able to address this issue more conclu- 
sively. It may also be interesting to note that in the re- 
lated bosonic problem, when a systematic expansion near 
four dimension is readily available, there are no charged 
critical points in the theory to the leading order 4^ This 
may also be contrasted with the well-known example of 
(albeit Lorentz invariant) scalar Higgs electrodynamics, 
for which the critical points, when they exist, are always 
charged j^^j^ 
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Probably the central message of this work is that pro- 
vided Lorentz invariance becomes emergent near critical- 
ity, for the Pz-orbitals well localized on carbon atoms the 
Lagrangian may be taken to contain only two (or with the 
physical spin, four) coupling constants. If there arc no in- 
tervening first-order transitions, one can infer then that 
there are two possible continuous metal-insulator tran- 
sitions, both governed by the same Gross-Neveu model 
in 2-1-1 dimensions, into the states that break either 
time-reversal or chiral symmetry. The residue of the 
quasiparticle pole on the metallic side plays the role of 
the metal's order parameter, and it vanishes continu- 
ously with a small critical exponent proportional to the 
fermion's anomalous dimension. In contrast, the specific 
heat coefficient 



lim — TT, 



(94) 



being dependent on the Fermi velocity only, at the tran- 
sition vanishes discontinuously from a finite value on the 
metallic side. Near the critical point and for the temper- 
atures much below the bandwidth we may assume that 
the specific heat obeys the scaling relation 



T 



(95) 



For small arguments the universal scaling function R(x) 
behaves as R{x) ~ 

2,2(^-1)/^^ SO that in the metallic phase 
one finds the usual quadratic temperature dependence 



(96) 



Recognizing the proportionality of the specific heat coef- 
ficient as (vp(ty)~'^ gives us yet another way to deduce 
Eq. (87). At the criticality, on the other hand. 



i?(oo). 



(97) 



with R{oo) expected to be finite. When z = 1, the spe- 
cific heat coefficient near criticality jumps therefore from 
i?(0) in the metallic phase to R{oo) at the critical point, 
and finally to zero in the insulating phascj^ 

Similarly, the optical conductivity near the metal- 
insulator transition will obey the scaling relation for t > 

with H{x) as a universal function, and with 

H{0) = J, (99) 

as the familiar universal dc conductivity per Dirac field 
in the metallic phase. In contrast to the specific heat, 
there is no (non-universal) dimensionful quantity such as 
vp in the scaling expression for conductivity, and conse- 
quently cr(0) is constant and universal in the entire metal- 
lic phase. Right at the transition then 



so the dc conductivity, while still universal, at the criti- 
cality should be different than in the metallic phase. Fi- 
nally, in the insulator the dc conductivity vanishes, so 
that the dc conductivity, similar to the specific heat co- 
efficient, in principle should show two universal disconti- 
nuities at the metal-insulator transition. 

One obstacle to experimental observation of these pre- 
dictions is that it has not been possible yet to tune the 
parameter ~(interaction/bandwidth) in graphene, and 
sample different phases of the system. The application 
of the magnetic field, however, changes this, since the ki- 
netic energy becomes completely quenched, and infinites- 
imal interaction immediately induces a finite gap. If the 
parameters of the system place it not too far from the 
metal-insulator transition, the gap m would obeyi^ 



(101) 



where I is the magnetic length, a = 1/A is the lattice 
constant, and t is the tuning parameter. G{x) is a (uni- 
versal) scaling function. The computation of the scal- 
ing function in the large-N limit, and the consequences 
of this scaling relation for experiment were discussed at 
length beforCfi^ Here we only wish to underline that the 
emergent Lorentz invariance of the metal-insulator criti- 
cal point, via its consequence that z = 1, implies precise 
proportionality between the interaction gap and the Lan- 
dau level separation at the criticality. 



m = vfVbG{Q), 



(102) 



(j{uj) — H{co)- 



(100) 



where G{0) is a universal number. Such a square-root 
magnetic field dependence of the gap is well known to 
arise from the long-range tail of the Coulomb interac- 
tions, but the above derivation serves to show that its 
origin may in principle lie in purely short-range interac- 
tions as well. 



X. SUMMARY 

We have presented the theory of electrons interacting 
via short-range interactions on honeycomb lattice, and 
in particular, determined the number and types of in- 
dependent quartic terms in the low-energy Lagrangian. 
Metal-insulator quantum critical points and the concomi- 
tant quantum critical behavior were discussed, with the 
particular attention paid to the consequences of the emer- 
gent Lorentz invariance. The minimal internally consis- 
tent local Lagrangian for spinless fermions is shown to 
contain only two Gross-Ncveu-like quartic terms. Gen- 
eralizations that would include long-range Coulomb in- 
teraction or spin of electrons were briefly considered. We 
also discussed the critical behavior of several key physical 
quantities on the metallic side of the transition, such as 
the Fermi velocity, the residue of the quasiparticle pole, 
specific heat, and the frequency dependent conductivity. 
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APPENDIX A: FIERZ IDENTITY 

For completeness, we provide the derivation of the 
Fierz identity. Assume a basis {T°',a = 1, ..16} in the 
space of four-dimensional matrices, and choose (r°)^ = 
pa _ (|pa)-i Then any Hermitean matrix M can be 
written as 



M = i(TrAfr")r'', 



(Al) 



with the summation over repeated indices assumed. This 
can be rewritten as 



and therefore it follows that 

1 



(A2) 



(A3) 



Applying this identity to the product of two matrix ele- 
ments then yields 



M,,iN, 



(A4) 



Finally, this leads to the expansion of a quartic term as 

(*(a;)M*(a;))(*(y)A^*(y)) = (A5) 

_ 1 (TrMF^iVF^ ) (i- (a;)r'' VI/ (y ) ) (^r (y )r- vfr (a;) ) 
16 

which is used in the text for x — y. The minus sign in 
the last line derives from the Grassmann nature of the 
fermionic fields. 



APPENDIX B: DIADIC FORM OF THE 
ASYMMETRIC MATRIX 

Any real A'^-dimensional matrix M can obviously be 
written as 



N 



(Bl) 



i=l 



where = {Mu, M2^, ...MNi), and (e,)j = S,j. In 
Dirac notation, 



N 



and 



N 



(B2) 



(B3) 



is the transposed matrix. There exists such a repre- 
sentation of the matrix M in any basis of vectors je^), 
as can be seen by multiplying M from the right with 

Let us now form a related symmetric matrix Sm = 
Ad^ M . Being symmetric, it can be written in the usual 
spectral form 



N 



Sm = y^^^J^i\n^){^x^ 



(B4) 



where {fii\^j) — Sij. We can now write, however, the 
matrix M in the particular eigenbasis of the associated 
symmetric matrix22. Sm 



N 



M = Y,\K^){fi^ 



(B5) 



From the definition of S and its spectral form we see that 
{Ki\Kj) = HiSij, and therefore 



AT 



(B6) 



where \Ki) — y^jz^i), and 



= 6ij. For a general 



asymmetric matrix, the basis /x and v are different, and 
the last equation generalizes the more familiar form for 
a symmetric matrix, where they are the same. 



APPENDIX C: RENORMALIZATION GROUP 
UNDER FIERZ CONSTRAINTS 

Here we provide an alternative formulation of the 
renormalization group transformation in presence of con- 
straints imposed by the Fierz identity. Let us demon- 
strate this method on the simplest example of the max- 
imally symmetric theory. Instead of choosing two in- 
dependent couplings and using Fierz transformation at 
intermediate stages of the calculation to transform any 
other generated quartic terms back into the chosen ones, 
one may use the kernel of the Fierz matrix to write the 
Lagrangian in terms only of the physical couplings from 
the outset, as in Eq. (53). The advantage of doing this is 
that no other quartic term besides the ones corresponding 
to the physical couplings can ever get generated then by 
the renormalization transformation. The set of couplings 
Al and A2 is therefore closed under renormalization. The 
computation to the quadratic order then yields 



dlnb 



d\nb 



-Xi- 24X1 -12X1X2, 



= -X2-72XI-AXI. 



(CI) 



(C2) 



The connection to Eqs. (63)-(64) in the text can be 
established as follows. Since the Fierz transformations in 
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this case imply 

= -3^2^ (C3) 

^Sl- 25*2, (C4) 

the Lagrangian Lint,max in Eq. (41) can obviously also 
be written as 

Lint,max = {gD2 - '^932 - "^gci) + {qai - gci) (C5) 

In the text we therefore have simply named the entire 
first bracket gD2, and the second gAi- But these can 



be recognized as particular linear combinations of the 
physical couplings Ai and A2: 

gD2 - 35S2 - 2gci - -2Ai - I2A2, (C6) 



(5Ai-5ci) = -2Ai. (C7) 

Such a connection is of course completely general, and in 
particular may be established between the three chosen 
couphngs in the Eq. (70) and the "physical couplings" 
determined by the vectors in Eqs. (58)-(60). 
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